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We argue that hadron multiplicities in central high energy nucleus-nucleus collisions are estab- 
lished very close to the phase boundary between hadronic and quark matter. In the hadronic picture 
this can be described by multi-particle collisions whose importance is strongly enhanced due to the 
high particle density in the phase transition region. As a consequence of the rapid fall-off of the 
multi-particle scattering rates the experimentally determined chemical freeze-out temperature is a 
good measure of the phase transition temperature. 
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The yield of (multi-) strange hadrons produced in cen- 
tral high energy nucleus- nucleus collisions was proposed 
two decades ago [lj as a signature of quark-gluon plasma 
(QGP) formation. Hadron yields observed in such colli- 
sions at AGS, SPS, and RHIC energies are found to be de- 
scribed with high precision within a hadro-chemical equi- 
librium approach [ISSIllESSHIIil; governed by 
a chemical freeze-out temperature T c h, baryo-chemical 
potential /i and the fireball volume \ c h- A recent review 
can be found in |ll| . Importantly, the data at SPS and 
RHIC energy comprise multi-strange hadrons including 
the Q and Q. Their yields agree with the chemical equi- 
librium calculation and are strongly enhanced as com- 
pared to observations in pp collisions. The time needed 
to achieve this equilibrium for Q baryons via two-body 
collisions was estimated pj to be much longer than rea- 
sonable lifetimes of the fireball. The observations were 
thus interpreted as a sign that the system had reached a 
partonic phase prior to hadron production [l^. flij . fl4| . 

In this note we argue that the chemical freeze-out tem- 
perature T c h is actually very close to the critical temper- 
ature T c of the QCD phase transition. This observation 
has a far reaching consequence: Since T c ^ has been mea- 
sured for different values of \x the approximate associa- 
tion of T c with T c /j implies that we have experimental 
knowledge of part of the critical line in the QCD-phase 
diagram. 

Let us first sketch our overall picture and detail our 
arguments subsequently. Hadro-chemical equilibration is 
achieved during or at the end of the phase transition. In 
particular, the number of strange quarks may be estab- 
lished in the plasma phase and/or hadronization of the 
QGP. During the very early stages of the hadronic phase 
the relative numbers of strange baryons and mesons 
K, K*, A, S, S, f2 are then realized at the thermal equi- 
librium values according to the Bose-Einstein or Fermi 
distribution 




/ p 2 dp{exp[(Ej(p) - fij)/T] ± (1) 
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Here E 2 — M? + p 2 and Mj proportional to the vac- 
uum mass of the hadron j, fij is the effective chemical 
potential, and gj counts degrees of freedom (for details 
see 0). The high accuracy of the distribution @J in re- 
producing the data suggests that T c h plays effectively the 
role of a universal temperature, governing simultaneously 
the chemical and kinetic distributions. 

In the hadronic picture the production of multi-strange 
hadrons can be described by multi-hadron strangeness 
exchange reactions 1 . The multi-hadron scattering is 
substantial, however, only in the immediate vicinity of 
the critical temperature T c . As T decreases the multi- 
particle rates drop very rapidly with a high power of the 
particle density. Below some temperature T c ^ very close 
to T c only two-particle interactions and decays remain 
as relevant processes for a change in the relative parti- 
cle numbers. These are too slow in order to equilibrate 
the distributions or to catch up with the decreasing tem- 
perature - chemical freeze-out occurs for T c h w T c . We 
proceed to discuss the three main points of this scenario 
in the following in more detail. 

(A) The QCD phase transition corresponds to a change 
in the effective degrees of freedom (from hadrons to 
quarks and gluons) in a narrow temperature interval 2 . 
For both the hadronic and quark-gluon phases sufficiently 
far away from T c the dominant processes in thermal equi- 
librium are two-particle scattering and decays. This is 



1 Production of multi-strange baryons by multi-particle collisions 
has also been considered by Il5l . Their argument focuses on 
anti-hyperon production at high baryon density, where indeed 
relatively short equilibration times are obtained. The authors 
conclude that their approach should not be applicable for RHIC 
energies, unless the hadronic phase has a rather long lifetime. 
Furthermore this approach does not take account of the expected 
rapid change of density near a phase transition which is central 
for our argument. 

2 Our central argument will make no distinction between a true 
phase transition and a rapid " crossover" . 
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FIG. 1: Comparison of the temperature dependence of energy 
density in lattice calculations with 2 light and 1 heavier quark 
flavor IT^I with results from the hadronic gas model of 0, 
with and without repulsive interactions. 

consistent with an effective (pseudo-)particle description. 
Close to T c , however, collective phenomena play an im- 
portant role. (As an example, near T c the c-resonance 
may behave like a particle with mass almost degenerate 
with the pions.) If T c h is close to T c the chemical equili- 
bration may be described in several equivalent pictures: 
multi-hadron scattering, time evolving classical fields or 
hadronization. We emphasize, however, that no picture 
should contradict a hadronic description. In the end, the 
chemical equilibrium distribution has to be established 
by fast processes involving hadrons (not quarks and glu- 
ons). 

Let us approach the phase transition (or a rapid 
crossover) from the hadronic phase. For T sufficiently 
below T c not much happens on the level of microscopic 
scattering processes between hadrons. The main effect of 
an increase of the temperature is an increase of the den- 
sity. Near T c , however, the density is so high that new 
dynamics can be associated with collective excitations. 
It is precisely the behavior of these collective excitations 
that triggers the transition. On the level of individual 
hadrons the propagation and scattering of collective ex- 
citations is expressed in the form of multi-hadron scat- 
tering. Since the collective dynamics becomes dominant 
only near T c the same holds for the multi-hadron scat- 
tering. 

We next argue that the temperature range where 
multi-hadron processes can dominate is actually very nar- 
row (typically a few MeV). This is due to (i) a rapid 
increase of the particle densities as a function of T and 
(ii) a very steep dependence of multi-particle scattering 
rates on the density of the incoming particles. Evidence 



for rapid energy density changes near T c comes from re- 
cent lattice QCD results This is shown in Fig. ^ 
where the observed rapid rise in energy density beyond 
the simple T 4 dependence reflects the large increase in 
degrees of freedom. The statistical model of a hadron 
resonance gas 0, exhibits only a moderate increase 
beyond the T 4 behavior (see Fig. determined by an 
interplay between the relevant number of hadronic states 
and the increased importance of repulsive interactions 
modeled by an excluded volume correction 3 . 

The rate of an individual multi-particle scattering pro- 
cess with riin incoming particles and average density n 
grows as n(T)™ iT1 . In addition, also the phase space 
increases with temperature. In the temperature re- 
gion of multi-particle scattering dominance one expects 
that many processes (with different rii n ) become of sim- 
ilar strength, thus increasing again the total rate. For 
the purpose of demonstration we model the tempera- 
ture dependence of a typical rate (say for scatterings 
which change the number of fi-baryons) by r ~ n(T) 7 , 
with 7 3> 1. Within a narrow temperature interval 
AT = 5 MeV this rate changes by a substantial factor 
(1 + with f3 = dlnn/dlnT > 1. This is demon- 

strated in Fig. |21 At T C) the multi-particle scattering 
rates are substantial. For T < T c , however, they drop so 
rapidly that only a small temperature interval around T c 
is left where the multi-particle scattering processes can 
dominate. The evaluation of these multi-particle rates is 
detailed in point (C) below. 

(B) Let us next turn to the issue of chemical equilib- 
rium and argue that two-particle scattering is insufficient 
to achieve or sustain it. We focus in the following on the 
analysis of data at RHIC energies and comment at the 
end on the applicability of our considerations at lower 
energies. The high accuracy of the statistical model pre- 
dictions in reproducing experimental particle ratios vary- 
ing over several orders of magnitude 0, EI should be 
taken as an indication that hadro-chemical equilibration 
has occurred in the system. Further, chemical equilibra- 
tion according to eq. Q requires two crucial ingredients: 
(i) The particle number changing reaction rates must be 
sufficiently high such that the particle distributions can 
adapt to a given T. (ii) At freeze-out the masses of the 
different hadrons must be proportional to their vacuum 
masses. From (ii) we conclude immediately that chemical 
equilibration and freeze-out occur in the hadronic phase, 
T c ft < T c . The particle distribution in the quark-gluon 
plasma has no memory of the individual hadronic vac- 
uum masses - for example, the relative density of strange 
particles (at the same given \i and T) would be deter- 



Neglecting the repulsive interactions the energy density of a 
hadronic gas diverges, reflecting the Hagedorn temperature (see 
dotted line in Fig. [I] and IH l 
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FIG. 2: Time tq = nn/rn needed to bring Q baryons into 
chemical equilibrium via multi-particle collisions. We display 
the dependence on the pion density and on T. The tem- 
perature scale is obtained as discussed in the text. The ar- 
row indicates the chemical freeze-out temperature T c h = 176 
MeV. 



mined by the strange quark mass m s rather than by the 
individual masses of K, K*,A, E, H, f2. 

To demonstrate that T c h is close to T c we exclude 
two possible alternative scenarios with T c ^ substantially 
smaller than T c : First, we show that an extended period 
of "hadronic chemical equilibrium" with T c h <T<T C is 
quite unlikely. Second, we generalize this argument to 
demonstrate that "late equilibration" at a temperature 
close to T c h but significantly smaller than T c is also dis- 
favored. 

The first argument is conceptually simple since the 
condition for chemical equilibrium at a temperature T 
substantially smaller than T c can be formulated by using 
the equilibrium rates and distributions. For our purpose 
complete thermalization of all quantities is not necessary 
- the " prethermalization" [18j of some rough quantities 
like relative particle densities and approximate momen- 
tum distributions is sufficient, such that it makes sense 
to speak about temperature and chemical potential (in 
the sense of eq. i[T]l) and to compute rates in thermal 
equilibrium. 

We first need the relevant time scale to which the two- 
particle scattering rates have to be compared. In equi- 
librium this is given by the inverse rate of decrease of 
temperature ry. For this purpose we assume, guided by 



recent results from two-pion correlation measurements 
two possible scenarios for the evolution 
of the fireball. In both cases we use the observation that 
the density decreases by only 30% between chemical and 
thermal freeze-out and our knowledge of T c h = 176 MeV. 
From the two-pion correlation data we obtain, for a cen- 
tral rapidity slice, the transverse and longitudinal radii 
at thermal freeze-out of 5.75 and 7.0 fm, a longitudinal 
expansion velocity /3|| =1, and a transverse expansion 
velocity f3±=0.5. The thermal freeze-out radii are to be 
understood as widths of Gaussian distributions and give 
a volume of V/= 3650 fm 3 . Scenario (1) assumes that 
the shape of the density distributions is the same (i.e. 
Gaussian) at thermal and chemical freeze-out and that 
accordingly an increase in density by 30 % is equivalent 
to a 30 % decrease in volume corresponding to V c h=2600 
fm 3 . In scenario (2) we use an initial volume of 1450 fm 3 
corresponding to a flat distribution over one unit of ra- 
pidity. With the assumption of isentropic expansion at 
the above velocities we find the longitudinal and trans- 
verse radius parameters at T c h as well as the duration 
of the expansion and the thermal freeze-out temperature 
Tf. For scenarios (1) and (2) the duration of the ex- 
pansion in the hadronic phase is Tf = 0.9 and 2.3 fm 
and the thermal freeze-out temperature is Tf = 158 and 
132 MeV. Consequently, in the hadronic phase near T c h 
the rate of decrease in temperature may be estimated as 
|t/T| = = (13 ± l)%/fm= (7.7 ± 0.6ml)- 1 . Note 
that these time scales are entirely consistent with the du- 
ration of pion emission, also obtained from the two-pion 
correlation data, of less than 2 fm. In both scenarios 
the lifetime of the fireball is rather short, leaving little 
room for an extended period of "hadronic cooking" at 
temperatures significantly below T c . 

As a simple example, a decrease in temperature by 
A T=5 MeV reduces the equilibrium ratio of il-baryons 
over kaons, n^/riK, by a factor F^k = 1.13 (using ther- 
mal model densities of Q). Adaptation of the particle 
distribution to the changing temperature requires Fqk — 
exp(aAi) = exp(arx AT/T ch ) with a = dlri(nn/nK)/dt. 
Let us define the rate of change of individual particle 
number densities as 



V 



+ n j V/V. 



(2) 



(The last term accounts for the decrease in particle num- 
ber densities due to the volume change.) Maintaining 
chemical equilibrium needs 



n K 
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AT 



(1.10-0.55)/fm = 0.55/fm. 

(3) 

The numerical evaluation of the two terms in the differ- 
ence agrees well with the direct evaluation of the term 
involving Fqk in eq. ©. For the difference in the rate 
of relative density change of 17-baryons to protons we ob- 
tain similarly a value of (1.10-0.90)/fm=0.20/fm. Both 
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rates are evaluated here at T c ^. Eq. Q can be obeyed 
either by the destruction of fi-baryons or the produc- 
tion of kaons. Since Sl-baryons decay weakly the effect 
of decays is completely negligible. A typical two-particle 
scattering ft + K — ► S + it (or + it — ► 3 + K) yields at 
most ffi/nQ = n^{va) =0.018 /frn. For this estimate we 
used a (large) strangeness exchange cross section of (7=10 
mb and a relative velocity of v—Q.Q. Similarly, for a typ- 
ical kaon production process it + it — > K + K one would 
obtain TkI^k = 0.18/fm, using the measured cross sec- 
tion of 3 mb j^. Clearly both numbers are much too 
small to maintain equilibrium close to T c h- Analogous 
arguments can be made for the f2/p ratio. The reaction 
7r + + £~ — > K~+p contributes to the rate of relative den- 
sity change of protons a value of 0.018/fm. One would 
need about 50 reactions with similar cross section (12 mb, 
[24|]) to keep the proton density in equilibrium. For the 
f2 baryons there is obviously no way to achieve this. We 
conclude that two-particle reactions and decays are not 
fast enough to maintain chemical equilibrium for multi- 
strange baryons in the hadronic phase near T c />. This 
finding is also supporte d by studies using cascade codes 
pi) and rate equations |26| . 

We next turn to the production rates for situations 
where thermal equilibrium has not yet been realized, 
e.g. during hadronization. The production rates should 
be consistent between different pictures for such pro- 
cesses and permit a hadronic description, at least towards 
the end of the chemical equilibration process. There- 
fore our picture does not require a detailed understand- 
ing of hadronization. Although our estimates of fj/rij 
have used thermal distributions for the incoming parti- 
cles our arguments can be extended to non-thermal situ- 
ations: there is no reason why the rates should be much 
larger or the available times longer. In the final approach 
to chemical equilibrium (needed in order to achieve the 
high accuracy of the observed thermal description of the 
data) the densities of incoming particles must already 
be close to equilibrium. Also only very rough features 
of the momentum distribution of the incoming particles 
are needed for an estimate of the magnitude of rj/rij. 
Therefore, the two-particle scattering is also too slow to 
achieve chemical equilibration in the production of mul- 
tistrange hadrons. In particular, this applies to a pos- 
sible picture of chemical equilibration by hadronization: 
hadronization at T much below T c would not lead to 
equilibrium abundances since there is not enough time to 
produce the multistrange baryons with rates correspond- 
ing to T.( Otherwise the hadronization picture would be 
in clear discrepancy with an equivalent hadronic picture 
for which the rates are dominated by two particle scat- 
tering.) This closes our argument: At T c h either multi- 
particle interactions must be important or the cross sec- 
tion must be dramatically larger than in the vacuum. 
Both possibilities are conceivable only for T c /, very close 
to T c . 



(C) The chemical equilibration of hadrons should be 
accessible to a hadronic description, at least in a rough 
sense. Consistency of the hadronic picture for equilibra- 
tion requires that multi-hadron processes changing the 
numbers of fl, £ etc. must be fast enough in order to 
build up the observed particle numbers at T c h- (We do 
not assume in this part thermal equilibrium with detailed 
balance of individual rates.) Keeping in mind the consid- 
erable quantitative uncertainties we now proceed to eval- 
uate rates for scattering processes involving more than 
two incoming particles and demonstrate the importance 
of such processes near T c . For an understanding of multi- 
particle interactions we write for the rate of scattering 
events per volume with m n ingoing and n out outgoing 
particles 

r(n m ,n out ) = fi(T) ni "\M\ 2 <j> (4) 

with 

*=n(/<^)<*> v (?4 (5) 

The rate is proportional to n;„ powers of the particle den- 
sities of the incoming particles that we denote for short 
by n(T). For the outgoing particles <fi is the Lorentz in- 
variant phase space factor which we evaluate case-by-case 
with the program given by |27^ . and which needs to be 
weighted (see below) with the thermal probability to find 
a particular cm energy in the initial state. The magni- 
tude of the squared transition amplitude 4 \A4 2 \ is evalu- 
ated using measured cross sections. We assume constant 
\M 2 \ independent of temperature and density. 

Inspection of measured cross section systematics shows 
that production cross sections of strange particles are 
at most a few mb, and usually much smaller, as is re- 
flected in the known strangeness suppression factor in 
hadronic collisions. Strangeness exchange reactions may 
reach cross sections in the 20 mb range. As cases in point 
we evaluate, using eq. Q, in the following the rate tq for 
f2 production through the reaction 2ir + 3K — > Nfl and 
similar rates for 2 and A production. 

For this case, 

m = nl(n K /n w ) 3 \M\ 2 <p. (6) 

The densities used to evaluate this expression are taken 
from the thermal model predictions which describe the 
measured yields at RHIC |t| using the excluded volume 



4 We use a normalization adapted to incoming fermions and out- 
going bosons. For each outgoing fermion the additional factor 
2E in the phase space integrals is absorbed here in the squared 
amplitude. For each incoming boson there is an additional factor 
2M in \A4\ 2 which is essentially cancelled by an additional factor 
l/(2Ei) inr. 
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correction. At T=176 MeV these are: = 0.174/fm 3 
and riK /n^—Q.172. Leaving out the excluded volume cor- 
rection would increase the densities by approximately a 
factor of 2. Note that n n and tik stand here for "generic 
pion and kaon densities" . Indeed, the incoming parti- 
cles can include all sorts of resonances like p etc. There- 
fore and tlk are not the thermal value of individual 
pion and kaon densities but rather comprise effectively all 
non-strange and strange degrees of freedom. The 2n'iK 
reaction is likely to be an important channel for SI pro- 
duction as it involves particles with the largest densities 
(pions, kaons) and the required amount of strangeness. 
Similarly we obtain EE's from 3tt2K and A's from AttK 
reactions. 

The numerical rate evaluation needs as input both the 
matrix element and phase space factor. The phase fac- 
tor 4> depends on yfs and needs to be weighted by the 
probability /(s) that the five-meson-scattering occurs at 
a given value of \fs. For this purpose we assume ther- 
mal momentum distributions for the kaons and pions in 
the entrance channel. The function / is evaluated nu- 
merically by a Monte-Carlo program. Its results were 
cross-checked for massless particles against an analytic 
evaluation 28]. The phase space factor is then obtained 
by folding with / in the energy range from threshold 
to infinity. 

For an estimate of the matrix element we note that 
the cross section for p + p — > 5ir has been measured. 
Close to threshold it takes a value of about 40 mb and 
is falling exponentially with cm energy -Js according 
to ctstt = 871mb • exp (- *fs 1.95/GeV) Q. For the 
SI + N —> 2?r + 3K reaction the threshold is 2.61 GeV 
and the peak of f(s) occurs at 3.25 GeV. To evaluate 
the corresponding matrix element we assume that the 
cross section at ^/s = 3.25 GeV is 6.4 mb, close to the 
pp — > 57r cross section at the same energy above thresh- 
old. From the known cross section and phase space the 
matrix element \A4\ 2 can be extracted by the usual for- 
mula. Using our normalization convention this yields 
\M\ 2 = 9.5 • 10 9 /GeV 8 . The final result for fi produc- 
tion through this channel is then tq = 1-39 • 10~ 4 /(fm 4 ) 
at T = 176 MeV and scales with the 5th power of the 
pion density. Furthermore, / scales approximately [2^ 
as -y/s 23 ^ 2 exp — y/s/T, leading to a further increase of tq 
with temperature (and hence density). 

We have alternatively evaluated the rate for SI produc- 
tion in a semi-classical approach, in which the standard 
two-body rate equation r = ninzicv) is generalized to 
multi-particle collisions. Inspired by the approach taken 
in cascade codes a reaction takes place if rii n particles 
approach within a volume V =47r/3-y/ (cr ine i Jn) . For the 
inelastic cross section we take a typical value of 0j ne ; = 40 
mb. A particular exit channel (such as SIN) is obtained 
by multiplying with a probability P^. For the reaction 
under consideration this yields 



r n = 8P x nUnK/n„) 3 yJ(*£/*3)(v). (7) 

Using for typical relative velocities (v) = 0.6 and P x = 
0.10 yields tq = 1.4 • 10~ 4 /fm 4 , indicating that a similar 
result as above is obtained with reasonable parameters. 

The meaning of this result is as follows: for a density 
n-x = 0.174/fm 3 , as used above, the final SI density of 
3.0T0~ 4 /fm 3 can be built up within a characteristic time 
tq = uq/tq = 2.2 fm. Since the SI yield scales in this 
approach as n 5 , and taking into account the temperature 
dependence of /, already an increase of n T to 0.2/fm 3 re- 
duces this time by more than a factor of three. The time 
tq is depicted in Fig. [3 as function of the pion density. 
We have also added a temperature scale in this figure. 
This scale is obtained as follows: we take the variation of 
energy density with temperature from the lattice results 
(see Fig.^l and assume that the density scales 5 as e/T. 
Fixing the overall scale by n n = 0.174/fm 3 at T = 176 
MeV as above this determines the T-dependence of n„ 
near T c . In consequence, close to T c , the SI equilibration 
time scales approximately as tq cx T -60 . 

We note that the above time of 2.2 fm is a reasonable 
time for hadronization and phase change, considering the 
necessary decrease by about a factor of 3 in the degrees 
of freedom and the concomitant volume increase. A de- 
crease of the pion density by 1/3 only will increase this 
time to about 27 fm; this time is even much longer than 
the total lifetime of the fireball as measured by two-pion 
correlations |l9| while we should consider here only the 
time between begin of hadronization and chemical freeze- 
out. 

We have checked the numerics of our approach by not- 
ing that, at equilibrium, detailed balance can be used to 
evaluate the 2ir3K — > Sin rate in the reverse direction 6 . 
Since both sides scale very differently with (overall) den- 
sity, setting both rates equal determines the equilibrium 
density (prior to resonance decays) for each temperature. 
The so-determined equilibrium density is to within 25% 
equal to that computed independently with our thermal 
model, lending strong support to our calculations. 

Our main results are fairly insensitive to the details 
of the calculation. As can be seen by inspecting Fig. 
the energy (and consequently particle) density is very 
rapidly increasing with temperature near T=176 MeV 
due to the phase transition. Density values of 20 % dif- 
ferent from those used in our calculations are reached 



5 Part of the increase in e is due to the increase in the number of 
effective degrees of freedom. Similarly, more relevant channels 
for five particle scattering open up. As argued above n v should 
be interpreted as an average weighted density n of all particles 
contributing to five particle scattering with Q-production. 

6 We thank C. Greiner and I. Shovkowy for pointing this out. 
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through temperature changes of less than 3 MeV. Pos- 
sible uncertainties in our rate estimates of even a factor 
of 2.5 would be compensated by such an increase or de- 
crease in density, indicating the stability of our estimate 
near T c . Furthermore, our rate estimates above are more 
likely overestimates because of the use of a comparatively 
large O N — * 3K2tt cross section. The resulting larger 
densities needed for chemical equilibration are, however, 
easily reached as in our picture the phase transition is 
passed through from above. 

We have, exactly along the lines for SI production, eval- 
uated the equilibration times for 5 and A production with 
values of = 0.71fm/c and ta = 0.66fm/c, indicating 
that all strange baryons have similar equilibration times 
with similar density and temperature dependencies. The 
corresponding time for protons and antiprotons is typi- 
cally shorter. We conclude that for all particle species 
the multiparticle rates are sufficient to produce the equi- 
librium abundances. 

It is an interesting question if for some particle species 
like the pion-proton-antiproton system the hadronic mul- 
tiparticle rates could be sufficient to maintain (restricted) 
chemical equilibrium even for some temperature range 
below T c . In this case the measured chemical temper- 
ature for the proton to antiproton ratio should reflect 
a lower T<T C ^, which is not the case observationally. 
This may be taken as an experimental indication that an 
extended period of hadronic equilibrium for T<T C is not 
realized. It may favour the idea that the relevant prether- 
malization process could be associated to hadronization. 

We note that thermal models have also been used 
[jfij to describe hadron production in e + e~ and hadron- 
hadron collisions, leading to temperature parameters 
close to 170 MeV. Indeed, this suggests that hadroniza- 
tion itself can be seen as a prethermalization process. 
However, to account for the strangeness undersaturation 
in such collisions, multi-strange baryons can only be re- 
produced by introducing a strangeness suppression fac- 
tor of about 0.5, leading to a factor of 8 suppression of ft 
baryons. In contrast to heavy ion collisions, tq exceeds 
here the available time. In the hadronic picture this is 
due to the "absence" of sufficient multi-particle scatter- 
ing since the system is not in a high density phase due 
to a phase transition. 

Has the critical temperature of the QCD-phase tran- 
sition been fixed by observation? The answer to this 
question needs a quantitative estimate of the difference 
AT = T c — T c h- An accurate determination is difficult 
since it involves the detailed understanding of equilibra- 
tion/ chemical freeze-out. From Fig. |2Jwe conclude that 
a temperature decrease of AT = 5 MeV below the criti- 
cal temperature lowers the five meson contribution to 
by more than a factor of 10. This factor is even larger if 
scattering processes with more than five incoming mesons 
dominate at T c h- Unless strangeness exchanging rates 
are unplausibly high at T c , such a sharp drop should 



make sufficient Q production impossible and we conclude 
AT-5 MeV. 

The accuracy of the experimental determination of T c 
could be limited by a possible temperature dependence 
of the hadron masses. Indeed, the hadronic yields deter- 
mine the ratio T c h/M rather than the absolute value of 
T c h- A large uncertainty in M(T) would reflect in a large 
uncertainty in T c h and we want to limit the size of this ef- 
fect. Using chiral symmetry arguments one may surmise 
that the masses of hadrons except for pions and kaons 
are proportional to the value of the order parameter a 
responsible for chiral symmetry breaking. As a result, 
they depend on T and /x even in the hadronic phase |3fll | 

M j (T) = h j (T,n)a(T,fj,) (8) 

Neglecting the T-and /^-dependence of the dimension- 
less couplings hj these masses scale proportional to 
cr(T, [i). However, not all hadron masses scale propor- 
tional cr(T, (i) - the pions and kaons scale differently. 
Therefore, too large mass changes are not easily consis- 
tent with the universality of chemical freeze-out. We no- 
tice that the "observed temperature" T b s = 176 MeV is 
fitted to the vacuum masses. The true freeze-out tem- 
perature T c h therefore obeys 

o-(T ch ,/x) = g(0,0) , g . 

Tch T Q bs 

Since a is expected to decrease for T> and \fi\ > (see 
e.g. [3(J ) we infer that the true freeze-out temperature 
T c /i ~ T c is lower than T h s . Already a moderate relative 
change of a by 10 % lowers T ch by 18 MeV [H| Keep- 
ing the rather vague character of our "error estimate" in 
mind we infer the critical temperature for small /i 

T c = 176^ 8 MeV. (10) 

A continuous crossover, second order phase transition 
or extremely weak first order phase transition make a 
ratio cr(T ch )/cr(0) close to one rather unlikely. The ob- 
served thermal distribution of the hadron yields strongly 
suggests, however, that the temperature dependence of 
the effective hadron masses should not be too substantial. 
This even could be interpreted as an experimental indica- 
tion that the QCD-phase transition is of first order! An 
accurate determination of the temperature dependence of 
hadron masses M(T)/M(0) by lattice simulations could 
greatly reduce the (lower) error in our estimate of T c . 

In summary, hadronic many body collisions near T c 
can concistently account for chemical equilibration at 
RHIC energies and lead to T c /j=T c to within an accu- 
racy of a few MeV. Any hadronic equilibrium scenario 
with T c h substantially smaller than T c would require that 
either multi-particle interactions dominate even much be- 
low T c or that the two-particle cross sections are bigger 
than in the vacuum by a large factor. Both of the lat- 
ter hypotheses seem unlikely in view of the rapid density 
decrease. 



7 



It was proposed 0, Q, ^ that the observed hadron 
abundances arise from a direct production of strange 
(and non-strange) particles by hadronization. How this 
happens microscopically is unclear. Nevertheless, in or- 
der to escape our argument that T c /j=T c one would have 
to argue that no (even rough) hadronic picture for this 
process exists at all - this is unlikely since the abundances 
are determined by hadronic properties (masses) with high 
precision. Second, one may question if the "chemical 
temperature" extracted from the abundances is a uni- 
versal temperature which also governs the local kinetic 
aspects and can be associated with the critical temper- 
ature of a phase transition in equilibrium. Indeed, in a 
prethermalization process, different equilibrium proper- 
ties are realized at different time scales. Nevertheless, all 
experience shows that kinetic equilibration occurs before 
chemical equilibration. It seems quite hard to imagine 
that chemical equilibrium abundances are realized at a 
time when the rough features of kinematic distributions 
(like relation between particle density and average ki- 
netic energy per degree of freedom) are not yet close to 
their equilibrium values. Defining the kinetic tempera- 
ture by the average kinetic energy one expects that the 
chemical and kinetic temperatures coincide at chemical 
freeze-out (with a typical precision at the percent level). 
More specific properties, like detailed balance of hadronic 
processes, may be realized later or never. 

The critical temperature determined from RHIC for 
T c /j T c coincides well with lattice estimates llfj for 
fj, = 0. The same arguments as discussed here for RHIC 
energy also hold for SPS energies: it is likely that also 
there the phase transition drives the particle densities 
and ensures chemical equilibration. The values of T c h 
and /i collected in would thus trace out the phase 
boundary for these energies. Whether the phase transi- 
tion also plays a role a lower beam energies is currently 
an open question. 

We thank J. Knoll for helpful discussions concerning 
phase space integrals. PBM and JS acknowledge the hos- 
pitality of the INT Seattle, where part of this work was 
performed. 
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